Quasiparticle localization in superconductors with spin-orbit scattering 
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We develop a theory of quasiparticle localization in superconductors in situations without spin 
rotation invariance. We discuss the existence, and properties of superconducting phases with local- 
ized/delocalized quasiparticle excitations in such systems in various dimensionalities. Implications 
for a variety of experimental systems, and to the properties of random Ising models in two dimen- 
sions, are briefly discussed. 



I. INTRODUCTION 

A powerful probe of the properties of a superconductor 
is obtained by studying the low temperature dynamics of 
the quasiparticles. In this context, we proposed [l| that 
all ground state phases of disordered superconductors can 
be characterized, at zero temperature, by their quasipar- 
ticle transport properties. The two general possibilities 
are that the quasiparticle excitations may be delocalized, 
analogous to a metal, or be localized analogous to an in- 
sulator. Previous papers P-ra have developed a theory 
of localization of quasiparticles in superconductors in sit- 
uations with spin rotation invariance. In this paper, we 
consider the case where the spin is not conserved. This 
may happen, for instance, in a singlet superconductor in 
the presence of spin-orbit scattering. Another example is 
provided by a triplet superconductor where the quasipar- 
ticles can exchange spin with the condensate, and hence, 
do not have conserved spin. Indeed, a number of su- 
perconducting systems, such as for instance the heavy 
fermion superconductors, are both strongly disordered 
and have strong spin-orbit scattering (and perhaps even 
triplet pairing). Thus, in order to understand the pos- 
sibility of quasiparticle localization in such systems, it 
is necessary to develop a theory that includes spin-orbit 
scattering in an essential way. Besides, by analogy with 
what happens in normal metals, spin-orbit scattering is 
expected to have profound effects on localization phe- 
nomena. 

Because the quasiparticle charge density is also not 
conserved in a superconductor, the only conserved quan- 
tity carried by the quasiparticles (at low energies) is the 
energy density itself. In the presence of impurity scatter- 
ing, the quasiparticle charge and spin densities in such a 
superconductor do not diffuse, as they are not conserved. 
Energy diffusion is possible, though. The corresponding 
transport quantity is the quasiparticle thermal conduc- 
tivity. 

From a theoretical point of view, quasiparticles in a 
superconductor with non-conserved spin are more appro- 
priately thought of as real (Majorana) fermions. Thus 
the problem we consider here is one of localization of 
Majorana fermions. While localization issues of com- 
plex (conventional) fermions have been explored in con- 



siderable detail, surprisingly, there has been very little 
theoretical work on corresponding issues for Majorana 
fermions. As we argue, the superconductor with non- 
conserved spin provides a natural experimental realiza- 
tion of such a system. We examine the possible phases 
(as characterized by quasiparticle transport) and the as- 
sociated phase transitions. It is of interest to distinguish 
between situations with and without time reversal sym- 
metry T, and we will consider each separately. (In the 
notation of Ref . Q , these correspond to Class Dili and 
Class D respectively). 

While the superconductor is our primary motivation, 
we note also that Majorana fermions arise in other con- 
texts as well. A well-known example is the two dimen- 
sional Ising model. There have been several studies of 
the properties of the two dimensional Ising model in the 
presence of randomness in the bond strengths, though 
there are still several poorly understood issues. The im- 
plications of this work for that problem will be considered 
briefly towards the end of the paper. 

We first show the existence, in two dimensions, of sta- 
ble "metallic" and "insulating" phases inside the super- 
conductor with delocalized and localized quasiparticle 
excitations respectively. The stability of the "metallic" 
phase in two dimensions has already been alluded to in 
Ref. M. We emphasize that both phases are supercon- 
ducting - they are distinguished by the nature of thermal 
transport due to the quasiparticles. These two phases are 
separated by a phase transition which is a novel "metal- 
insulator" transition inside a superconductor. In Ref. [y, 
we discussed possible experimental realizations of such 
phase transitions. The universal critical properties of 
this transition depend, of course, on whether time rever- 
sal is a good symmetry or not. In T-non-invariant sys- 
tems, the insulating phases may be further characterized 
in terms of their values of the Hall thermal conductance. 



The dimensionless ratio 






is quantized in units of 



1/2 . Phases with different values of this quantized Hall 
thermal conductance are topologically distinct, and are 
separated by phase transitions. In Figure H, we show a 
schematic phase diagram. Note that, in the case with no 
T, apart from the "metal-insulator" transition, there are 
also transitions between the "metal" and the "quantum 
Hall" phase, and from the insulator with Kxy/T = to 



the "quantum Hall" phase. Further, there is a multicriti- 
cal point, where all three phases come together. In three 
dimensional systems, "metallic" and "insulating" phases 
are again possible with the transition between the two 
being in a new universality class for localization. In the 
rest of the paper, we will substantiate these claims, and 
analyse the properties of each phase in further detail. 
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FIG. 1. Schematic zero temperature phase diagram for the 
two dimensional superconductor in the absence of both spin 
rotation and time reversal invariances. 



II. MODELS AND GENERAL FORMALISM 

Consider a general lattice Hamiltonian for the quasi- 
particles in a superconductor with strong spin-orbit scat- 
tering: 



n = j:^cl. 
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(1) 



Here i,j refer to the sites of a lattice, and a, (3 are spin 
indices. We assume that t and A are both short-ranged 
in space. 

We will focus on two cases - with and without time 
reversal symmetry (T). If present, time reversal is im- 
posed through an antiunitary "time-reversal" operator T 
which transforms the c operators as: 



r 









(2) 



Note that T<t*T^^ = cr (with a a vector of Pauli ma- 
trices) so that the electron spin is odd under time re- 
versal. With time reversal invariance present we require 
that 



TTiT-^ = n. 



(3) 



Note that H has no special symmetries (other than 
possibly T). In particular, neither the charge nor spin is 
conserved. It is convenient then to work with Majorana 
fermions rjua , rj2ia defined through 



-^ {Vila +irji2a) 



(4) 



The Hamiltonian when expressed in terms of the f] 
Fermions takes the form 



n 



v'^Ht], 



(5) 
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with 77 = rjiaa, (a, b — 1,2) and H = H'^"''"''" is a matrix 
in (ij), (a/?), ab space. By definition, H is hermitian and 
moreover H = —H, so that H is pure imaginary. Thus 
the problem of quasiparticle localization in a supercon- 
ductor with spin-orbit scattering is, in essence, one of 
localization of Majorana fermions. In particular, in the 
case where even T is not a good symmetry, the Hamil- 
tonian Eqn. is the most general one describing non- 
interacting Majorana fermions in a disordered system. 

Time reversal symmetry is easily imposed on the Ma- 
jorana Hamiltonian Eqn. 0. Under the action of the 
antiunitary operator T, it is readily seen that 



TViaa = i{cFy)al3{T^)abriibl3- 



(6) 



Here a,b = 1,2 and r is a Pauli matrix in ab space. The 
time reversal invariance condition Eqn. ^ combined with 
Eqn. |6| implies 



<7yTzH*CFyT2 = H. 



(7) 



We note in passing that for a superconductor of spin 
polarized Fermions (appropriate, say, in the A~ 1 phase 
of superfluid 3 — He) in the absence of time reveral in- 
variance, the general Hamiltonian can still be written 
in the form Ti = rf^ H-q, with H an antisymmetric and 
pure imaginary matrix in position (i,j) and "particle- 
hole" {ah) space. But for spinlcss quasiparticles the con- 
dition for time reversal invariance is different, leading 
to different symmetries. Specifically, with time rever- 
sal invariance the lattice Hamiltonian depends on a real 
symmetric hopping matrix, t^, and a real antisymmet- 
ric gap matrix, A^j. In this case, when re-expressed in 
terms of Majorana fermions, the Hamiltonian becomes 
Ti = rf" Hrj with H = try + i/^Tx- This Hamiltonian ma- 
trix can equivalently be expressed as H = iAr^ + c.c, 



with 



iTy and A — t 



A an arbitrary real 



matrix. This "off-diagonal" form is very different than 
the two cases with spinful electrons, and in fact belongs 
to a different symmetry class - a class studied by Gade 
and Wegner S|. Henceforth, we focus exclusively on the 
spinful case. 

We are interested in understanding the nature of en- 
ergy transport (in both spinful cases, with and without 
T) by the excitations described by the Hamiltonian. For 
this purpose, it is actually convenient to adopt the follow- 
ing trick. We consider two identical copies of the system. 
To describe the two copies, we introduce two Majorana 
fields ri and C, and consider 



H = rfHri + CHC- 



(8) 



It is now possible to combine the 77 and C, fields into a 
single complex fermion /: 



1 



(9) 

(10) 

(11) 



Note that the number of / particles is conserved. It is 
then possible to consider transport of this conserved /- 
number density. This may be quantified by a conductiv- 
ity af. As the Hamiltonian Ti. describes non-interacting 
/-particles, the thermal conductivity Kf of the / particles 
is related to u/ by a Weiedemann- Franz law (as T ^ 0): 
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(12) 



Since the Hamiltonian Ti. represents just two identical 
copies of the original system described by H, it is clear 
that the thermal conductivity (k) of the rj particles is 
exactly half that of the /-particles: 



«/ 



(13) 



Combining this with Eqn. n3, we see that calculation of 
K is reduced to calculation of af . Similarly, all the ther- 
modynamic properties of Ti. may be obtained by halving 
the corresponding property calculated with Ti. 
It is convenient to define the /-Green's function 



G 
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(14) 



where e is a positive infinitesimal and E is the en- 
ergy (measured from the Fermi energy). The condition 
H* = —H immediately implies that 



G*{E) = ^G{-E), 



(15) 



where we have suprressed all the indices on G. In partic- 
ular, for states at the Fermi energy, we have 



G*{E^O) = -G{E^O), 



(16) 



so that G{E = 0) is pure imaginary. This Green's func- 
tion can be convenintly expressed, as usual, as a func- 
tional integral over Grasmann variables /ma, /ma with 
the action 



S = ifiH - ie)f. 



(17) 



The symmetry Eqn. Hq implies that this same generating 
functional can be used to calculate transport quantities 



III. SYSTEMS WITH TIME REVERSAL 
SYMMETRY 



The discussion has so far been completely general. We 
now specialize to the case with T symmetry. 



A. "Metallic" Phase 



We first consider the situation in which the disorder is 
weak. We assume that there is a finite, non-zero mean 
free path le set by the impurity strength in such a manner 
that the /-particle motion is diffusive on larger scales. (In 
terms of the original physical system, this corresponds to 
diffusion of energy). We are interested in describing the 
effects of quantum interference on this diffusive motion. 
We assume also that, in the absence of quantum interfer- 
ence effects, the density of states at the Fermi energy of 
this diffusive system is finite and non-zero. In that case, 
it is possible to follow standard techniques to derive a 
replica non-linear sigma model field theory to describe 
the physics at length scales large compared to the mean 
free path Ig. As the procedure is sufficiently well-known, 
we merely state the results. The field theory is described 
by the action 



S= f d'^x^Tr (VO)^(VO)) - eTr{0 + O^ 



), (18) 



where 0{x) £ 0{2n) is a 2n x 2n orthogonal matrix- 
valued field with n the number of replicas. When e = 
the action has a global 0{2n) x 0{2n) symmetry, O — > 
AOB with A and B orthogonal matrices, which is broken 
down to the diagonal 0{2n) [A^ = B) by the e term. The 
coupling constant g is inversely proportional to the bare 
/-conductivity ct9 (i.e the conductivity on the scale of the 
mean free path). The density of quasiparticle states at 
the Fermi energy is given by 



p=|^(rr(0 + 0^)). 



(19) 



Consider a renormalization group transformation 
where short distance fluctuations are integrated out, and 
the coordinate x is rescaled as x — *■ a;' = xe~^ . The lead- 
ing quantum interference corrections to diffusion can now 
be obtained from the known p| perturbative (3 function 
of this field theory in the replica limit. The result, in two 
dimensions, is 



dg_ 

dl 



47r' 



(20) 



Note that g decreases as I is increased. Thus the per- 
turbation theory (in powers of g) gets better at large 



length scales. For a system of size L, at zero tempera- 
ture, we may integrate the flow equation upto a scale I* 
given by leC^ ~ L to get 



9iL) 



50 
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(21) 



where go is the bare value oi g. For large L, this therefore 
gives 

An 



9{L) 



in'u 



(22) 



Thus g{L) goes to zero logarithmically with the system 
size. As Uf is inversely proportional to g, it follows that 
Uf diverges logarithmically with the system size. At finite 
temperature, in an infinite system, it is natural to expect 
that the quantum interference effects will be cut-off at a 
finite dephasing length scale L^ ~ T~'p due to interaction 
effects not included in the model. We therefore have 



'^f 



ln( - 
T 



(23) 



at the lowest temperatures. As this also determines the 
thermal conductivity k(T) of the original system, we have 



K 

f 



In 



(24) 



The considerations above establish the existence of a 
"metallic" phase with delocalized quasiparticle excita- 
tions in two dimensions in the model Hamiltonian Eqn. 
H describing the superconductor in the presence of time 
reversal invariance, but no spin rotation invariance. In 
striking contrast to normal metals, quantum interference 
effects also lead to singular corrections to the density of 
states in a superconductor Q . In our previous work Q , 
we demonstrated this in the spin-rotation invariant cases. 
We now show that the density of states is enhanced in 
the situation considered in this paper. In particular, we 
show that in two-dimensions it actually diverges in the 
thermodynamic limit. 

To see this, consider the action Eqn. O at finite e. 
The density of states is obtained from 



P_ 
Po 
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(25) 



£=0 



where J^ is the free energy density defined through 
exp(-L''j^) ^ I VO exp(-5). 



(26) 



The flow of e under the rcnormalization group may be ob- 
tained straightforwardly. To leading order in g the result 



ll 



(2+ A) 

V 8ttJ 



(27) 



This is readily integrated to give 



e(0 = e(0)exph/ + i-y gil')dA 



(28) 



The free energy density scales according to 

n9(0),e(0)) = e-2'^(5(0,e(0), (29) 

where we have specialized to two-dimensions. For a 
system of size L, run the RG till a scale I* such that 



L 



L. The density of states is then 



p{L) = - 



-21* 



An 



dT{g{r),e{r)) 
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After scaling out to I* the mean free path is comparable 



to the system size, so that — ^ — 3e(i*) 
Therefore, 



p(L)^exp(^^y 



dl'gil') 

For large Z, g(V) « ^. Thus, we have 
p(L)- Win 



const.. 



(31) 



L 



(32) 



Thus the density of states at the Fermi energy diverges. 
The behaviour of the density of states as a function of 
energy in an infinite system may also be found similarly 
(or simply guessed from the equation above) to be 



P{E) 



■4- 



(33) 



The Einstein relation can now be used to infer that the 
heat diffusion constant D also diverges on approaching 
the Fermi energy as 



D{E) 



In- 



1 

'e' 



(34) 



Thus the "metallic" phase has an infinite heat diffusion 
constant at zero temperature in two dimensions. 

The divergence of the density of states has obvious con- 
sequences for the low temperature thermodynamic prop- 
erties of this phase. For instance, the specific heat be- 
haves as 



C(T)^T\/lni 



(35) 



In three dimensional systems, the stabihty of the 
"metalhc" phase can be estabhshed by simple power- 
counting arguments. Quantum interference effects are 
then irrelevant at long length scales. The thermal con- 
ductivity then goes to zero linearly with the temperature 
at low T: 



H3D 

T 



const. 



(36) 



The density of states at the Fermi energy is finite, and 
non-zero. However, quantum interference effects do lead 
to a singular y'lE] cusp in the density of states as a func- 
tion of the energy (See Ref. M for an analogous discussion 
in the superconductor with conserved spin) so that 



p{E) - p{E - 0) 



(37) 



Note that the density of states increases with decreasing 
energy as also happens in d = 2. 



B. "Insulating" phase 

For strong disorder, it is possible to have a phase with 
localized quasiparticle excitations in any finite dimen- 
sion. In the terminology of Ref. ||I|, this is a supercon- 
ducting " thermal insulator" . In this phase, the ratio y 
goes to zero with the temperature. The density of quasi- 
particle states also goes to zero at the Fermi energy. To 
see this, consider the lattice Hamiltonian Eqn. [ij in the 
extreme limit where the / particles are localized to a sin- 
gle site. The Hamiltonian for a single site is constrained 
to be of the form H — auyTx + bry with a, b real. This 

Consider now 



has two eigenvalues given by ±-\/a^ 4- 
the case where the distribution of a, b has finite, non-zero 
weight at a = 6 = 0. Then, the density of states p{E) 
averaged over the disorder vanishes as |-E|. If the dis- 
tribution has vanishing weight at a = 6 = 0, then p{E) 
vanishes faster than linearly. Including hopping between 
the sites should not change this result so long as we are in 
the localized phase. (As with the superconductors with 
conserved spin, having a finite density of states requires 
a diverging weight at a = 5 = which is presumably 
unphysical, and non-generic). We thus conclude that the 
density of quasiparticle states vanishes, at least as fast as 
\E\ in the localized phase. 



IV. TIME REVERSAL BROKEN SYSTEMS 

We now move on to systems without time reversal 
symmetry. As mentioned earlier, this corresponds to 
studying the general localization properties of a non- 
interacting disordered system of Majorana fermions. It 
is therefore appropriate to call the phase with extended 
states a "Majorana metal" , and the phase with localized 
states a "Majorana insulator" . 



A. Majorana metal 



To address issues such as the stability and properties 
of the Majorana metal, it is useful to think in terms of a 
replica sigma model field theory which permits a system- 
atic study of quantum interference corrections to diffusive 
energy transport. This is readily done using standard 
techniques. The result is 



S 



d?x—Tr{dQf 
25 



eTr {(TyQ) 



(38) 



where Q — V^ CyV is a 2n x 2n matrix, V is an 0(2n) 
matrix and n is again the replica index. When e = 0, the 
action is invariant under the group 0{2n), Q -^ W^QW, 
with W an orthogonal matrix. With non-zero e, invari- 
ance of the action requires W = exp{iG) with G pure 
imaginary and antisymmetric of the form, W = Say+iA. 
Here 5* is an n by n real symmetric matrix and A is n by 
n real antisymmetric. Since S + iA is n hy n Hermitian, 
the group has evidently been broken down to U(n) by 
the e term. Thus S describes a non-linear sigma model 
on the manifold 0{2n)/U{n). The coupling constant g 
is again inversely proportional to the (longitudinal) /- 
particle conductance. The density of states of the quasi- 
particles may be obtained from 



1 
2^ 



(rr(a,Q)) 



(39) 



Actually, on symmetry grounds, there is a topological 
term allowed in the action. This follows from the observa- 
tion that H2 (0(2n)/C/(7i)) = Z is non-trivial. This will 
be important to understand the behaviour of the thermal 
Hall conductivity later. 

The perturbative RG flows of the couplings g and e in 
two dimensions in the replica limit are 



(40) 
(41) 



We may use these to extract the physical properties of 
the Majorana metal phase in exactly the same manner 
as in the previous section. We therefore simply state the 
results. First note that the (marginal) irrelevance of the 
coupling g implies the stability of the Majorana metal in 
d — 2. As the temperature goes to zero, the (longitudi- 
nal) thermal conductivity k^x in a finite system of size L 
behaves as 



dg 
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(42) 



In an infinite system at finite T, where tfie quantum inter- 
ference is cut-off by dephasing due to interaction effects, 
we liave 



K. 



XX T / -L 

— ~ in — 
T \T 



(43) 



Similarly, the density of states at the Fermi energy now 
diverges logarithmically with the system size: 

p(L)~lnf. (44) 

I'e 

As a function of energy in an infinite system, we have 



p(S)~lni 



(45) 



This leads, for instance, to a specific heat which depends 
on temperature as 



cm 

T 



In- 
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(46) 



Since Hxx/T and the density of states diverge in the same 
manner with energy, determinining the behavior of the 
thermal diffusion coefficient would require going to sec- 
ond order in the perturbative RG. The diffusion coeffi- 
cient could then also diverge, but not more rapidly than 
double logarithmically with energy. 

In three dimensions, the behaviour of the supercon- 
ducting "metal" phase is qualitatively the same irresepec- 
tive of whether time reversal symmetry is present or not. 
Therefore, the discussion in the previous section of the 
three dimensional case applies here as well. 



C. Major ana quantum Hall phase 

In the absence of time reversal symmetry, there is po- 
tentially another transport property which can be used 
to distinguish zero temperature phases: the thermal Hall 
conductivity Kxy In particular, in the Majorana insula- 
tor in two dimensions, the ratio -^ approaches quantized 
values as the temperature goes to zero. Phases with dif- 
ferent values of the quantized thermal Hall conductance 
are topologically distinct and are separated by phase tran- 
sitions. In recent work with X.G. Wen [Q, we studied 
the physics of this novel quantum Hall system in some 
detail. 

From the point of view of the replica non-linear sigma 
model field theory discussed earlier in this section, the ex- 
istence of insulating phases with quantized thermal Hall 
conductance can be attributed to the presence of a topo- 
logical term in the action. This takes the form 



Si 



top 



'x-—Tr (QdxQdyQ) . 

IDTT 



(47) 



The quantity multiplying 9 is equal to i x integer for any 
given configuration of the Q field. Thus the partition 
function is periodic under 9—^9 + 2TTm for any inte- 
ger m. Physically, just as in the Pruisken field theory 
for the conventional integer quantum Hall transition, 9 
is proportional to the (bare) value of the /-particle Hall 
conductivity a^ . As shown in Ref. jl^l, this is in turn 
proportional to the ratio of the thermal Hall conductivity 
to the temperature. 



V. PHASE DIAGRAM AND TRANSITIONS 



B. Majorana insulator 

At strong disorder, in any dimension, it is possible to 
find phases where the quasiparticle excitations are local- 
ized. The (longitudinal) quasiparticle heat conductivity 
i^xx goes to zero rapidly with the temperature, in such a 
phase. We will call this the Majorana insulator, as this 
corresponds to a localized phase of Majorana fermions. 
The density of states of the Majorana insulator may be 
found by an argument analogous to the one in Section 
[II B. In brief, we consider a single site Hamiltonian for 
two Majorana species (corresponding to a single spinless 
complex Fermion) which is constrained to be of the form 
H = aTy with a real and random. If the distribution of a 
has a finite, non-zero weight at a = 0, then there is a finite 
density of states at zero energy. We again expect this to 
hold through out the localized phase, so that the density 
of states at the Fermi energy can be non-zero. Note that, 
like in a conventional insulator in non-superconducting 
systems, a density of states that vanishes at the Fermi 
energy is also possible. 



We now discuss the (zero temperature) phase diagram 
and phase transitions for the dirty superconductor with 
non-conserved spin. Again, we consider the cases with 
and without T separately. 



A. Time reversal invariant systems 

As argued in the previous section, "metallic" and "in- 
sulating" phases are possible in both d = 2 and d = 3. A 
schematic phase diagram as a function of the bare ther- 
mal conductivity is as shown below in Figure 0. 

SC- THERMAL METAL SC-THERMAL INSULATOR 



V. 



FIG. 2. Schematic zero temperature phase diagram for the 
three dimensional superconductor in the presence of spin orbit 
scattering. 



The critical points separating the two phases are in new 
universahty classes for Anderson localization. There are 
some interesting differences between the scaling proper- 
ties of the two and three dimensional systems. In three 
dimensions, the ratio y , the thermal diffusion coefficient, 
D, and the density of states po at the Fermi energy are 
finite constants in the "metal" and are zero in the "insu- 
lator" . It is therefore natural to expect that these will go 
to zero continuously (with some universal critical expo- 
nent) as the transition is approached from the "metallic" 
side. In two dimensions however, in the "metal" , these 
quantities are all infinite in the limit of zero temperature 
and infinite system size. (They are zero in the "insula- 
tor"). Exactly at the transition point, conventional scal- 
ing arguments ||ll| imply a constant value for y. The 
behaviour of the density of states and diffusion coeffi- 
cient is more unclear. It is perhaps natural to suggest 
that these will be a constant, though we certainly cannot 
rule out other possibilities. 



B. Time reversal non-invariant systems 

We first consider the three dimensional case. The 
phase diagram is similar to that above for the time rever- 
sal invariant system. However, as the density of states 
may be non-zero in the "insulator" , we expect that it 
generically is non-zero at the transition point as well. 

The two dimensional case has a richer phase diagram 
due to the possibility of quantum Hall phases. For sim- 
plicity, we show, in Figure nl only the two insulating 

phases with K^y = and -^ = i (in units of -jjf-)- 
Note that there are three distinct phase transitions, and 



a multicritical point in this phase diagram. First consider 
the Majorana metal-insulator transition. Both -^^ and 
Po are infinite in the Majorana metal. In the Majorana 
insulator, -^^ is zero while pq may be non-zero. As in 
the T invariant case, conventional scaling arguments im- 
ply a finite value for ^^ at the transition point. For the 
density of states, it is natural to guess that it is finite at 
the transition, though again it is hard to rule out other 
possibilities. 

Now consider the transition from the Majorana metal 
to the quantum Hall phase. This transition is in the 
same universality class as the transition from the Majo- 
rana metal to the insulator with Kxy — 0. The quantum 
Hall phase is a Majorana insulator with -^ = i. The 
distinction between these two insulating phases is due 
to the presence of (heat) current carrying edge states - 
these are expected to be unimportant in determining the 
properties of the transition to the metal. Below, we will 
provide a more formal argument in support of this claim. 

Finally, consider the transition between the insulator 
with Kxy = and the one with -^ = 5- ^ theory for 
this transition is obtained by considering a particular re- 



alization of the two phases. Consider a system of spinless 
fermions paired into a p^ + ipy superconducting state in 
two dimensions. In Ref. [|0[, it was shown that such 
a superconductor has a quantized thermal Hall conduc- 
tivity -^^ = i if the chemical potential is positive (the 
"BCS" limit), and has Kxy = if the chemical poten- 
tial is negative ("molecular" limit). Thus, as the chemi- 
cal potential is varied through zero, there is a transition 
from the Majorana quantum hall phase to the Majorana 
insulator. Ref. [in] also examined the theory for this 
transition, and argued that, at least at weak disorder, 
it is correctly described by a theory of relativistic Majo- 
rana fermions with random mass. It is well-known that 
the random mass is irrelevant at the pure free Majorana 
fixed point [[7| . Thus the critical theory is known exactly 
in this case. If we make the important assumption that 
there is a unique fixed point describing the Majorana 
insulator-quantum Hall transition, then the arguments 
above identify it with the free relativistic Majorana fixed 
point. There is, however, some reason to question this 
assumption - see Section VI. 



Some more insight into the phase diagram and the 
transitions comes from considering the properties of the 
replica non-linear sigma model field theory describing 
the quasiparticles in the two dimensional superconductor 
with spin orbit scattering, and no T. This is described 
by the action Eqn. pq supplemented with the topological 
term Eqn. B^. The small g regime is described by the 
perturbative calculations of the previous section. Note 
that the presence of the topological term does not af- 
fect those results. Indeed, the parameter 9 plays no role 
in perturbation theory, and does not renormalize. Thus 
there is a line of fixed points in the g, 9 plane at g = 
with 9 arbitrary. As the value of 9 is proportional to the 
value of ^, this ratio varies continuously in the Majo- 
rana metal. 

Now consider the behaviour of the sigma model at 
large g. It is natural to expect that, in this limit, the 
physics is captured by a strong-coupling expansion, and 
that the resulting phase corresponds, physically, to a lo- 
calized phase. (For conceptual purposes, it may be con- 
venient to think in terms of a lattice regularized version of 
the field theory.) First consider 9 = 0. The correspond- 
ing localized phase has Kxy = 0. Similarly, the insulator 
with 9 = 2% has ^^ := 2- Just as in other sigma model 
field theories with topological terms, the large g phase 
for < < TT is expected to be continuously connected 
to the large g phase at 9 = 0, i.e it is an "insulator" with 
Hxy — 0. Similarly the large g phase with n < 9 < 2tt 
is expected to be continuously connected to the large g 
phase with 9 = 27r, i.e an "insulator" with ^^ = i 
The transition between the two localized phases occurs 
ai 9 = TT. The small g "metallic" phase is separated 
from these localized phases at large g by a phase bound- 
ary. The symmetry 9 —> 2t: — 9 implies that this phase 
boundary be symmetric about the 9 = n line. 
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FIG. 3. Phase diagram of the sigma model field theory 
describing the two dimensional superconductor with no spin 
rotation or time reversal invariances. 



In the field theory, the transitions between the small g 
"metal" and the large g "insulator" at = 0, 27r are in the 
same universality class as the physics is invariant under 
9 -^ 9 + 2tt. What about the transition from the "metal" 
to the "insulator" at other values oi 9 7 We suggest that, 
for 9 ^ n, these are in the same universality class as the 
transition at 9 = 0. This then implies that Kxy is con- 
tinuous across the "metal-insulator" transition. Finally, 
the point where the "metal-insulator" phase boundary 
crosses the 9 = n line in the phase diagram of the field 
theory corresponds to the multicritical point where the 
"metal" and the two "insulating" phases come together. 



VI. RANDOM BOND ISING MODEL IN D = 2 



In this section, we briefly discuss the possible implica- 
tions of this paper to the finite temperature properties of 
classical random bond Ising models in two dimensions. It 
is well-known J7l that this system admits a description in 
terms of non-interacting disordered Majorana fermions. 
Specifically, the partition function for the Ising model can 
be written as a functional integral over Majorana fields 
that live in two dimensions. This, in turn, may be reinter- 
preted as a generating function for the wavefunctions of 
a quadratic quantum Hamiltonian of Majorana fermions 
in two spatial dimensions. This Majorana Hamiltonian 
has no special symmetries - the system is therefore in the 
universality class of the two dimensional superconductor 
with neither spin rotation nor time reversal invariance. 



It is therefore natural to try to identify various phases 
of the Ising model with corresponding (T = 0) phases in 
the two dimensional superconductor with no spin rota- 
tion or time reversal invariance. First consider the pure 
Ising model. Both the high and low temperature phases 
correspond to gapped phases in terms of the Majorana 
fermions. Thus these correspond to insulating phases 
of the Majorana fermions. Disorder in the Ising model, 
obtained by making some of the bonds random, would 
tend to fill up the gap. However, if the disorder is weak, 
the resulting low energy states will be strongly localized. 
Thus, at weak disorder, both the high and low temper- 
ature phases of the Ising model correspond to localized 
phases of the Majorana fermions. How then do we dis- 
tinguish between the two? The distinction is topological, 
with the quantized value of the thermal Hall conductance 
differing by 1/2 between the two phases. Strong support 
for this suggestion is obtained by examing the properties 
of the system near the transition between the two phases. 
In the pure Ising model, the critical point is described by 
a massless relativistic Majorana theory. Moving off criti- 
cality introduces a mass for the fermions with the sign of 
the mass distinguishing the two phases. The result that 
the dimensionless Kxy jumps by 1/2 at the transition can 
now be established by direct calculation. The critical 
theory is identical to the one describing the transition, 
at weak disorder, in the px + ipy superconductor of spin- 
less fermions between the BCS and molecular limits. We 
thus identify the Ising transition at weak disorder with 
the thermal quantum Hall transition. 

Now consider introducing strong (bond) randomness 
into the Ising model. If all the bonds are still positive, 
the transition is believed to be always controlled by the 
pure Ising fixed point. The situation is more interesting if 
some of the bonds are made negative. As the concentra- 
tion of negative bonds is increased, there is a multicritical 
point (the Nishimori point, see Figure^) ||l2[ on the phase 
boundary, beyond which the transition is no longer con- 
trolled by the pure Ising fixed point. The properties of 
this multicritical point have been the focus of a number 
of investigations [Q over the last many years. Despite 
this, there is no detailed understanding of the theory of 
this point, or of the properties at the phase boundary at 
lower temperatures. What does the the Nishimori point 
correspond to in the language of the Majorana fermions? 
There appear to be two distinct possibilities - we will 
discuss evidence supporting either below. 

(i) For the problem of localization of Majorana 
fermions in two dimensions, the assumption that there 
is a unique fixed point controlling the thermal quantum 
Hall transition leads to it's identification with the free 
relativistic Majorana theory. This assumption, which is 
perhaps natural in the fermion language, then points to 
the identification of the Nishimori point with the multi- 
critical point separating the Majorana metal, insulator, 
and the quantum Hall phase. If true, then it is natu- 



ral to expect that the ferromagnetic phase transition in 
the Ising model at low temperatures below the Nishiniori 
point is actually the Majorana insulator-metal transition. 
Some evidence in support of this is provided by the nu- 
merical results of Ref. [H which found no signs of local- 
ization of the fermions in the non-ferromagnetic phase at 
low temperature close to the phase boundary. On the 
other hand, note that the paramagnetic phase at high 
temperature and weak disorder corresponds to a local- 
ized phase. If the scenario outlined above is correct, then 
we are led to infer the existence of a finite temperature 
phase transition outside the ferromagnetic phase in the 
Ising model associated with a delocalization of the Ma- 
jorana fermions. It is unclear what this transition means 
in the Ising language, and even whether it happens at 
all. A natural candidate would have been a transition to 
a spin glass phase at low temperature - however there is 
strong numerical evidence for the absence of spin glass 
order at finite temperature in two dimensional Ising sys- 
tems. It seems possible that a delocalized ("metallic") 
phase of Majorana fermions would correspond to a phase 
in which both the Ising spin and the dual Ising disorder 
parameter are simultaneously zero. 
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FIG. 4. Phase diagram of the two dimensional Ising model 
with random bond strengths, p is the concentration of nega- 
tive bonds, and T is the temperature. The point A'^ is the 
Nishimori multicritical point. SG refers to the spin glass 
phase which is believed to exist only at zero temperature. 



(ii) A different scenario is obtained by dropping the 
assumption that there is a unique fixed point control- 
ling the thermal quantum Hall transition. Instead, for 



the Ising model, we assume that the only finite temper- 
ature phase transition is associated with the destruction 
of the ferromagnetic order. If the concentration of nega- 
tive bonds is large enough to destroy the ferromagnetism 
at zero temperature, we assume that the resultant state 
is a spin glass. These assumptions are perhaps reason- 
able expectations for the Ising model. Then, the finite 
temperature ferromagnetic and paramagnetic phases are, 
in the fermion language both localized phases which are 
distinguished by their thermal Hall conductance. The ex- 
istence of the Nishimori point then implies the existence 
of a multicritical point (at strong disorder) in the phase 
boundary between these two localized phases. Then, the 
thermal quantum Hall transition will be controlled by the 
free relativistic massless Majorana theory at weak disor- 
der, but by a different fixed point at strong disorder. 
The spin glass phase in the Ising model presumably cor- 
responds to the Majorana metal. However, it may seem 
a bit puzzling, in this scenario, why the Majorana metal 
which is a stable phase of disordered Majorana fermions 
in two dimensions is only realized in a set of measure zero 
(a line at T = 0) in the Ising phase diagram. 

Which one of these two scenarios is actually realized 
and the resultant consequences both for the Ising model 
and the Majorana localization problem, we leave as an 
intriguing open question. 



VII. DISCUSSION 

In this paper, we have discussed the physics of local- 
ization of quasiparticles in a superconductor in situations 
where the quasiparticle spin is not a good quantum num- 
ber. Our discussion was entirely based on models of non- 
interacting quasiparticles - including the effects of inter- 
actions is an important and interesting issue which we 
leave for future work. Even within the non-interacting 
theory, we have found a rather rich phase diagram, and 
a number of as yet unexplored phase transitions. Several 
experimental systems to which this work is of relevance 
can be imagined - here we briefiy mention a few. 

(i) A particularly attractive prospect for probing the 
physics discussed here is in a Type // superconductor 
in a strong magnetic field in the presence of spin-orbit 
scattering impurities. At low fields, the system will be 
in the superconducting "insulator" phase. With increas- 
ing field, under conditions discussed in Ref. ||l|, there 
will be a transition (at zero temperature) to a supercon- 
ducting "metal" phase. This transition can be probed, 
for instance, by measurements of the low temperature 
quasiparticle heat conductivity W. 

(ii) A number of recent experiments p5| have mea- 
sured low temperature heat transport by the quasiparti- 
cles in the heavy fermion superconductors. We note that 
as these systems typically have strong spin-orbit scat- 
tering, the results of this paper are of potential impor- 



tance. Some of the heat conductivity measurements have 
been motivated by the possibihty of identifying the cor- 
rect pairing symmetry in these superconductors. Such 
identification could be seriously hampered by the local- 
ization issues discussed here. In particular, if the disorder 
is strong enough to localize the quasiparticles, it is hard 
to infer whether the pure system has a gap to quasiparti- 
cle excitations or not from heat transport measurements 
alone. 

(iii) The properties of superfluid He — 3 in porous me- 
dia have been the subject of some experimental stud- 
ies |16 . In this context, it is interesting to ask if the 
fermionic He — 3 quasiparticles in the superfluid are lo- 
calized or delocalized at zero temperature. This may, 
perhaps, again be probed by heat transport experiments. 
It also interesting to consider the properties of superfluid 
He — 3 on a disordered two dimensional substrate. In the 
A phase, at weak disorder, a thermal quantum Hall effect 
|lO| is predicted. With increasing disorder, it is possible 
(though not necessary) that there is a zero temperature 
phase transition where the quantization of the thermal 
Hall conductivity is destroyed before the superfluidity is 
destroyed. If this happens, this would be an experimen- 
tal realization of the thermal quantum Hall transition 
discussed in the previous section |10| . 

Finally, the general problem of localization of Majo- 
rana fermions arises in other contexts as well - a specific 
example being the random bond Ising model in two di- 
mensions as discussed briefly in Section W\ Progress in 
the localization problem may therefore provide a route 
to further our understanding of the random bond Ising 
model. 
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